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INTRODUCTION 


Five  paper*  dealing  with  practical  aspect*  of  theoretical  work  done 
at  the  Applied  Physics  Laboratory  on  the  subject  of  the  aerodynamic  char¬ 
acteristics  of  wings  in  supersonic  flow  are  presented  in  this  BUMBLEBEE 
report.  These  papers,  which  hare  appeared  as  internal  memoranda  of  the 
Applied  Physics  Laboratory,  represent  the  combined  efforts  of  APL/JHU 
personnel  and  outside  contributors  as  noted  in  the  references. 

The  first  paper  derives  lift  coefficients  and  center  of  pressure 
location  for  fiat  plates  of  polygonal  planform  from  fundamental  consid¬ 
erations,  using  Busemann* s  conical  field  method.  The  other  papers  utilize 
Busemann* s  second  order  approximation  formula  to  determine  the  aerodynamic 
characteristics  of  certain  types  of  wings  having  finite  thickness.  Con¬ 
sequently  there  will  be  some  overlapping  of  results,  but  no  material  is 
deleted  inasmuch  as  the  two  methods  are  quite  different.  Attention  is 
called,  however,  to  differences  in  nomenclature  between  the  first  and 
the  remaining  papers.  A  short  discussion  of  the  sectional  properties 
of  various  airfoil  shapes  and  the  optimum  type  is  given  at  the  end  of 
this  report. 

It  is  believed  that  these  papers  cover  most  of  the  unclassified 
engineering  work  on  the  effect  of  wing  planforms  available  at  this  time. 
The  work  is  by  no  means  complete,  there  being  many  practical  planforms 
for  which  no  theory  has  been  developed,  particularly  in  the  case  of  wings 
with  finite  thickness.  Accordingly  the  table  on  the  following  page  ha* 
been  constructed  tc  indicate  the  types  for  which  information  can  be  found 
in  this  report,  to  show  where  it  can  be  found,  and  to  indicate  the  mist¬ 
ing  data.  In  this  table  the  Roman  numeral  refers  to  the  paper  in  this 
report  as  indexed  in  the  Table  of  Contents,  whereas  the  Arabic  numeral 
refers  to  the  page  on  which  the  expression  for  the  particular  coefficient 
may  be  found.  An  X  indicates  that  the  exact  information  is  not  yet  known. 
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Lift  i 8  two-dimensional  and/or  center  of  pressure  is  at  centroid  of  area. 


I.  APPLICATION  OP  BUSSMABJP  S  COHICAL  PIELD  METHOD  TO  THIS  WINGS* 


by  Robert  M.  Snow 


Introduction 


In  this  paper  the  pressure  distribution,  lift  coefficient  and  center 
of  pressure  are  determined  for  several  plane  wings  of  polygonal  planform 
at  small  angles  of  attack,  by  the  method  of  "conical  fields".  The  effect 
of  a  dihedral  bend  is  also  obtained.  Squares  and  products  of  perturbation 
velocities  have  been  neglected,  since  the  method  is  based  on  the  Prandtl- 
Glauert  linearised  potential  equation.  Modifications  due  to  wing  thick¬ 
ness,  viscous  effects,  and  interference  effect  (with  a  fuselage  or  with 
other  wings)  are  likewise  neglected. 

The  method  of  conical  fields  in  supersonic  aerodynamics  was  developed 
by  Busemann  (Ref.  1),  who  applied  it  to' several  important  problems.  Stew¬ 
art  (2)  has  solved  the  problem  of  a  delta  wing  by  essentially  the  same 
method.  A  conical  field  corresponds  to  a  linear  homogeneous  solution  of 
the  iinearised  (or  Prandtl-Glauert)  potential  equation  for  supersonic 
flow: 


A 


(Mz-I) 


=o 

az2 


a) 


Here  M  is  the  Mach  number  of  the  main  stream,  which  is  moving  along  the 
z  -  axis.  The  perturbation  velocity  components  (u,  v,  w)  are  also  solu¬ 
tions  of  Eq.  (l)  and  are  homogeneous  of  degree  sero,  i.e.,  u,  v,  and  w 
are  constant  along  any  ray  emanating  from  the  origin.  The  particular 
simplicity  of  conical  fields  lies  in  the  fact  that  after  a  transformation 
the  perturbation  velocity  components  (u,  v,  w)  are  obtained  as  solutions 
of  Laplace’s  equation  in  two  variables.  Busemanu  credits  this  transfor¬ 
mation  to  Chaplygin,  (3)  who  made  use  of  it  in  a  formally  similar  problem. 
Since  the  general  theory  has  buen  discussed  recently  by  Stewart  (2)  it  is 
only  needed  to.  state  the  principal  result  in  the  form  in  which  it  will  be 


*  This  paper  is  a  revision  of  CM-265  originally  published  as  an  internal 
memorandum  of  the  Applied  Physics  Laboratory. 
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Aerodynamic  Characteristics  of  Wings  at  Supersonic  Speeds 


utilized.  I<6 i  ,Ai  he  the  Mach  angle  (sin"^  l/M),  and  x,  y,  z  a  rectan¬ 
gular  coordinate  system  with  the  2  -  axis  pointing  downstream,  and 

I  "V  v2  a .  y  2 

®=  tarf  y/x  r  =  — A=tan 

The  transformation 


is  such  that  homogeneous  functions  of  degree  zero  which  satisfy  Eq.  (1) 
also  satisfy  Laplace’s  equation  in  the  polar  coordinates  r,  p  ). 

The  evaluation  of  the  streamwise  component  (w)  of  perturbation  velocity 
is  of  primary  importance  since  the  aerodynamic  forces  on  the  wing  are 
determined  by  w  alone.  This  follows  from  the  linearized  Bernoulli 
equation, 


p  =  p-^Wuj 


(3) 


which,  like  Eq.  (1J,  results  from  neglecting  squares  and  products  of 
perturbation  velocities  in  the  corresponding  exact  equation.  In  many 
problems  of  this  type,  including  those  considered  here,  w  may  be  deter¬ 
mined  without  further  reference  to  the  components  u  and  v. 


Application  of  Busemann*  8  Conical  Field  Method  to  Thin  Wings 
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Characteristic  Cones,  Boundary  Conditions 


To  illustrate  the  type  of  boundary  condition  needed  to  determine 
the  conical  field,  consider  the  special  case  of  a  rectangular  wing. 

This  rectangular  wing  may  he  regarded  as  the  result  of  cutting  off  the 
ends  of  a  two-dimensional  airfoil.  This  operation  causes  a  modification 
in  the  flow  (originally  two-dimensional),  which  may  he  referred  to  as 
the  "tip  effect".  In  this  connection,  a  fundamental  distinction  should 
he  made  between  subsonic  flow  and  supersonic  flow.  For  subsonic  flow 
(differential  equation  of  elliptic  type)  the  tip  effect  dies  off  asym- 
pototically  with  increasing  distance  inboard.  For  supersonic  flow  (dif¬ 
ferential  equation  of  hyperbolic  type)  the  tip  effect  falls  to  zero  at 
a  certain  finite  distance,  and  the  entire  effect  is  contained  within 
a  region  bounded  by  real  characteristic  surfaces.  For  linearized  super¬ 
sonic  flow,  the  domain  of  influence  of  any  point  is  bounded  by  a  "Mach 
cone",  which  is  one  nappe  of  a  cone  opening  downstream  with  semi -vertex 
angle  equal  to  the  Mach  angle /x  .  Figure  1  represents  a  section  by  a 


plane  perpendicular  to  the 
main  stream.  The  Mach  cones 
from  the  tips  of  the  leading 
edge  divide  this  plane  into 
three  types  of  regions.  In 
the  central  region  (I)  the 
flow  is  in  all  respects  the 
same  as  if  the  wing  were  of  . 
infinite  span,  because  no 
point  of  this  central  region 
lie 8  in  the  domain  of  in¬ 
fluence  of  any  point  removed 
in  the  mental  process  of  ob¬ 
taining  the  rectangular  wing 
from  an  airfoil  of  infinite 
span.  On  the  other  hand  the 
perturbation  velocity  compon¬ 
ents  are  zero  in  the  exterior 
region  (III),  because  no 
point  of  this  region  lies  in 
the  domain  of  influence  of 
any  point  on  the  rectangular 


m 


I  Flow  Some  at  for  about 
Infinite  Spon  Airfoil 

n  Conical  Field 
in  Free  Stream 


FIG.  I 

FLOW  FOR  RECTANGULAR  WING, 
SECTION  BY  PLANE  PERPENDICULAR 
TO  STREAM 


wing.  The  requirement  of  con¬ 
tinuity  leads  to  boundary  con¬ 
ditions  which  must  be  satisfied  by  the  perturbation  velocity  components 
u,  v  and  w  on  the  boundary  of  each  conical  transition  region  (II).  On 
the  boundary  between  regions  (II)  and  (III),  u,  v  and  w  must  vanish. 

Or.  the  boundary  between  regions  (I)  and  (II),  u,  v  and  w  take  on  the 
(constant,  two-dimensional)  values  of  region  (I). 


Similar  statements  apply  also  to  the  example  of  a  swept-back  lead¬ 
ing  edge  (Fig.  2).  The  two  lines  forming  the  leading  edge  are,  of  course, 
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finite  in  the  actual  case,  hut 
the  effect  of  their  finiteness 
will  he  confined  hy  character¬ 
istic  cones  passing  through 
the  points  at  which  the  leading 
edge  changes  direction.  The 
possibility  of  treating  more 
general  polygonal  wings  hy 
this  method  follows  from  these 
remarks,  for  the  polygons  that 
may  he  treated  hy  this  method, 
there  are  isolated  regions  of 
uniform  flow  (identical  with 
the  flow  for  an  infinite  wing 
at  a  certain  angle  of  yaw)  sep¬ 
arated  hy  regions  of  transition 
in  the  Mach  cones  which  start 
from  each  vertex  of  the  polygon 
and  open  downstream.  It  is 
convenient  in  the  following 
discussion  to  refer  to  these  • 
special  Mach  cones  simply  as 
"the  Mach  cones"  or  "the  Mach 
cone". 


I,  Flow  »ome  ot  lor  Infinite  ipof.  airfoil  at  onglt  of  yow  Tf/2  -  (, 
It  Flow  torn*  at  for  Infinite  tpon  oirfoll  ot  onglt  of  yew  TT/2-S, 
II  Coniccl  Field 
IH  Frit  Stroom 


FIG.  2 

FLOW  FOR 
SWEPT-BACX  LEADING  EDGE 

Boundary  conditions  must 
also  he  given  over  that  part  of 
the  wing  which  lies  inside  the 

Mach  cone.  The  velocity  component  normal  to  the  wing  must  have  the  same 
value  as  on  the  rest  of  the  wing,  because  of  the  condition  that  there 
he  no  flow  through  the  wing.  The  boundary  condition  for  w  is  that  the 
normal  derivative  of  w  he  sero  on  the  wing.  This  follows  from  the  re¬ 
quirement  of  irro.tationality  (implicit  in  the  use  of  the  potential  $  ), 
so  that  Sur  _  dv 

by  dz 

and  the  condition  that  v  is  constant  on  the  wing.  The  fact  that  the  wing 
does  not  lie  exactly  in  the  plane  y  -  0  is  neglected;  this  has  no  effect 
on  the  first  order  perturbation.  This  simplification  is  made  throughout, 
so  that  the  single  of  attack  enters  only  in  the  boundary  values  and  not  in 
the  position  of  the  boundaries. 


Rectangular  Wing 


The  problem  of  a  plane  rectangular  wing  was  successfully  solved 
hy  Busemann  (1),  hy  essentially  the  method  employed  here.  Schlichtlng 
(4)  had  previously  considered  the  same  problem  hy  a  difierent  method, 
and  obtained  a  false  result  because  of  an  analytical  error.  When  this 
error  is  corrected,  Schlichting’ s  method  becomes  consistent  with  the 
conical  field  method,  not  only  for  rectangular  wings,  hut  also  for  the 
raked  wings  considered  in  the  next  section.  (Baked  wings  were  not  treated 
by  pither  of  these  authors,  hut  their  methods  are  eaeily  extended  to  this 

case . ) 


Application  of  Busemann1 b  Conical  Field  Method  to  Thin  Wings 

Attention  may  "be  confined  to  one  end  of  the  wing.  The  origin  of 
coordinates  is  taken  at  the  end  of  the  leading  edge.  The  positive  x  - 
e±ie  is  directed  spanvise  away  from  the  wing.  The  y  -  axis  is  normal 
to  the  wing;  the  direction  of  the  positive  y  -  axis  may  he  regarded  as 
"upward".  The  positive  s  -  axis  is  in  the  direction  of  flow  (the  co¬ 
ordinates  being  such  that  the  wing  is  considered  to  he  at'  rest  with  the 
air  flowing  past  it).  The  velocity  of  the. incident  flow  is  W  ,  the 
angle  of  attackflC  ,  and  the  Mach  angla/U  «  sin~l  1/M.  For  a  wing  of 
infinite  span,  the  usual  tvo- dimensional  theory*  gives 

U7  =  ocW  tan  p.  - 

above  the  wing  and 

VJ  =  -  VJn 

below  the  wing.  In  terms  of  the  polar  variables  R  and  <p  ,  the  bound¬ 
ary  conditions  fprb'j'  are: 


bvj 

—  =  0 

for  ?=iTT.  0<R<A 

dtp 

*  # 

ijj  -  +W& 

for  75<?<Tr,  R  =  A 

uj«  0 

for  +g- ,  R=  A 

U*  “UJa, 

for  -tf  <  >  R  =  A 

These  boundary  conditions  of  course  apply  equally  well  in  the(^,<p) 
plane,  replacing  the  circle  E  -  A  by  the  circle  r  s  1.  To  find  w,  it 
is  only  necessary  to  solve  Laplace's  equation  in  the  polar  coordinates 
r  and<p  with  the  above  boundary  conditions.  The  solution  may  be  con¬ 
structed  by  the  usual  methods  from  the  particular  solutions  r8  sin  s<p. 
From  the  first  boundary  condition  it  is  seen  that  S  must  be  half  of  an 
odd  integer.  The  moat  general  function  which  is  harmonic  in  the  cut 
circle,  which  satisfies  the  required  conditions  along  the  cut,  and' which 
is  an  odd  function  of  <p  (as  required  by  the  symmetry  of  the  boundary 
conditions)  is 


sin(n+ifeH 


*  See,  for  example,  von  Misss  and  Friedrichs,  Fluid  Dynamics^  p.  237, 
Brown  Univercity,  1942. 
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This  is  not  the  conventional  Fourier  series  for  a  function  which  is 
periodic  with  period  2/r?  5  however,  it  is  useful  to  notice  that  the 
series  is  periodic  with  period  ^TT  .  This  suggests  the  extension  of 
the  problem  to  a  two-sheeted  Rienann  surface  by  analytic  continuation 
across  the  cut.  It  is  convenient  to  retain  the  symmetry  by  consider¬ 
ing  (p  to  range  from  -&*?T  to  +  JUT  ", 

then  the  boundary  values  on  the  arcs  in  the  "lower**  sheet  of  the  Riem&nn 
surface  (that  is,  -^fT<(p<-7r  '7T<  <p<  Z'lf) 


are  to  be  assigned  in  such  a  way  that  the  even  harmonics,  which  cannot 
appear  in  the  series  for  w  drop  out.  This  is  done  by  assigning  bound¬ 
ary  values  in  the  lower  sheet  which  are  the  negative  of  the  values  at 
corresponding  points  of  the  tapper  ("physical")  »beet.  Then 


(-l)n4i0a»  sin(fH-^a)‘^2 
TT  2n  +1 


Using  the  symbol  R  to  indicate  "real  part  of",  the  Fourier  series  for 
w  may  be  written  as 


UJ  = 


*n  +  /« 


rr 


n*o 


2n  +  I 


UJ  _  i  2V7  co$(*'2-\)  I  2V7cos(^2  A) 
—  -  ^tan  — - 7-37— - —  tan  1 - - 


UJ, 


l-r 

For  <P  =  =  tan 


rr 


uj  _  2  rzr 


l-r 


00 


TT 


l-r  Tt 


sin 


l  /j£L. 

V  I  +  rs 


(4) 


(5) 
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On  the  lower  surface  of  the  wing 

(9  =  -7t; 


w  is  merely  changed  in  sign.  The  right  member  of  Eq.  (5)  or  (51  )  is 
shown  as  a  function  of  R/A  in  Fig*  3*  From  the  form  (5')i  it  is  clear 
that  the  curve  is  symmetric  with  respect  to  the  point  (1/2,  1/2),  so 
that  the  average  value  of  w  along  any  spanwise  line  from  R  =  0  to  R  =  A 
is  just  Woo / 2.  The  triangle  cut  off  by  the  Mach  cone  contributes  just 
one  half  as  much  to  the  lift  as  an  equal  area  located  in, the  region  where 
two-dimensional  calculations  are  valid,  (it  will  be  seen  in  the  next 
section  that  the  same  simple  result  holds  for  raked  wings  with  leading 
and  trailing  edges  perpendicular  to  the  flow.)  Formulae  for  lift  co¬ 
efficient  and  center  of  pressure  are  now  easily  obtained.  If  Cj,  co 
is  the  lift  coefficient  for  a  plane  wing  of'  infinite  span,  and  the 
lift  coefficient  for  a  rectangular  wing  of  total  span  s  and  chord  c, 

=  ^  (6) 

L-oo 

The  center  of  pressure  is  located  at  a  distance  back  from  the  leading 
edge  expressed  by 

c  i-2/3C/S  tan  >JL  (7) 

l-l/2C/s  tanjx 


The  above  discussion  tacitly 
assumes  that  the  two  Mach  cones  from 
the  tips  do  not  intersect  on  the 
wing.  However,  this  restriction  is 
seen  to  be  unnecessary.  If  is 
the  (disturbance)  potential  inside 
one  of  the  cone's,  the  potential 
inside  the  other  cone,  and$co  the 
potential  for  a  wing  of  infinite 
span  with  leading  edge  perpendicu¬ 
lar  to  the  stream,  then  in  the 
region  common  to  the  two  cones  the 
potential  is  $  =  * -§<*,  .  To 

verify  this,  it  may  be  noted  that 
cik<{  $oo  are  solutions  of 
the  Prandtl-Glauert  Eq.  (1),  and 
since  that  equation  is  linear,  is 
a  solution  as  well.  Also  $  and  its 
first  derivatives  (u,  v,  w)  are  con¬ 
tinuous  across  the  conical  surfaces 
bounding  the  region  in  question. 

It  may  be  said  that  the  Htip  effects” 
from  the  two  tips  are  additive, 
since  the  equation  defining  may 
be  written 


A  is  the  froctionol  distonce  inboord,  spanwise, 
from  the  wing  tip  to  the  Mach  cone. 

FIG.  3 

PRESSURE  DISTRIBUTION  ALONG  SPAN 
FOR  RECTANGULAR  WING  OR  FOR 
RAKED -BACK  RECTANGULAR  WING 
WITH  161  <M 
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It  should  he  noted  that  the  flow  in  the  region  in  question  is  not  a  "conical 
field*,  The  combination  of  the  two  conical  fields  with  different  vertices 
i 8  never  a  conical  field,  although  it  approaches  a  conical  field  asymptoti¬ 
cally  downstream. 

This  solution  is  limited  by  the  condition  that  the  Mach  cone  from 
one  tip  should  not  intersect  the  other  end  of  the  wing.  If  this  happens 
a  further  alteration  is  needed  to  satisfy  boundary  conditions  at  the  edge, 
and  the  difficulty  of  the  problem  is  increased  enormously. 

Since  overlapping  "tip  effects*  are  additive,  the  lift  coefficient 
and  center  of  pressure  are  easily  calculable  as  long  as  the  Mach  cone  from 
one  wing  tip  does  not  cut  the  other  wing  tip.  After  calculation,  it  is 
found  that  Eqs.  (6)  and  (7)  also ‘hold  for  the  overlapping  case.  It  may 
be  noticed  that  as  the  quantity £ tan  u  increases  from  0  to  1  (the  highest 
value  for  which  the  solution  applies;,  the  ratio  C^/Cj,  ^  decreases  from 
1  to  1/2,  and  the  center  of  pressure  moves  from  C/2  forward  to  C/3. 


Baked  Wing 


The  following  treatment  applies  tc  wings  with  a  positive  rake  angle 
0  ,  defined  in  such  a  way  that  the  leading  edge  is  greater  than  the 
trailing  edge.  If  £>/•<  ,  the  problem  is  quite  simple.  The  pressure 
on  the  wing  is  uniform  and  is  to  be  calculated  from  two-dimensional 
theory.  This  is  seen  by  considering  the  wing  as  carved  out  of  an  infin¬ 
ite  span  wing,  which  only  involves  removing  portions  whose  domain  of 
influence  does  not  contain  any  part  of  the  final  trapezoidal  wing. 


if  o<S<>u  ,  the  boundary  conditions  differ  from  the  previously 
considered  special  case  of  the  rectangular  wing  only  in  that  the  "cut" 
or  radial  line  on  which  hi cr/gris  Q  (and  across  which  w  is  necessarily 
discontinuous;  does  not  extend  as  far  as  the  coordinate  origin.  The 
cut  along  <p*±7T  in  the  (R,  (p  )  plane  runs  from  R  =  A  =  tan ju  to  R  - 
+  tan  $  •  In  the  (r,  <p  )  plane  the  c  it  is  still  on  the  ray  (psdb'rT 
and  runs  from  r  3  1  to 


If  «  + 


ton  jjl 
ton  S 


i*  fo 


The  new  problem  may  be  solved  by  finding  a  conformal  transformation 
which  carries  it  back  to  the  already  solved  problem  of  the  rectangular 
wing.  To  do  this  it  is  advantageous  to  map  both  the  new  and  the  old 
problem  onto  a  half -plane.  Letting  ,  the  transformation 


£=f-Vpr 


or 


maps  the  upper  half  of  the  circle  m\  onto  the  upper  half  plane 
S  q  .  The  transformed  boundary  conditions  relate  to  the  real 
axis  £ 
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VJ  = 

VJ  - 

VJ  = 
OUJ 

bn  “ 


0 

0 

WoD 

0 


for  S,>0 
for  = 

for  - !  <  §,,<  0 
for  §,<§^<-1 


ton  m  _  ^ 
ton  $ 


The  ultimate  objective  is  to  find  a  transformation  which  connects  the 
£,  plane  with  the  plane  corresponding  to  a  fictitious  rectang¬ 

ular  wing.  The  £'  plane  is  mapped  onto  the  %l  plane  by  the  same  trans¬ 
formation  which  connects  £  and  §  .  The  boundary  conditions  in  the' 
plane  are  obtained  from  those  in  the  §  plane  by  setting  §o=  -  oo 

The  transformation  which  carries  the  £  plane  into  the  4'  plane  it  now  clear; 
it  must  be  1>he  homographic  transformation  which  carries  the  points 


into 


(to,-1'0) 


(oo,-!,  o) 

This  is  easily  found  to  be 


c/  (S»+t)t 

S.-S 


/ 

The  transformation  from  the  £  plane  to  the  £ 


p”?ane  is  now  known. 


since 
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By  eliminating  ^  and  §  it  in  possible  to  obtain  £  as  a  function 
of  £  .  In  the  plane  of  the  wing,  where  <p=0  or  7T  and  §  =  A /Rq*-1 
the  result  simplifies  to 


Rei<p'_  Re^-fRe 
A  ^  A-Ro 


0  or  TT 


Here  R0,  are  the  polar  coordinates  of  the 
0  or  TT  according  as  the  rake  is  positive  or 
of  the  wings 


R.'_  R-R. 
A  A-  R. 


edge  of  the  wing;  <p0  is 
negative.  On  the  surface 


As  R*  goes  from  0  to  A,  R  goes  from  R0  to  A,  and  the  dependence  of  R* 
on  R  is  linear.  Thus  the  lift  for  a  raked  wing  also  varies  according 

*°  fW-£  sin'1  VT 

II 


where  A  *  ft* /A  is  the  fraction  of  the  spanwise  distance  from  the  edge 
of  the  wing  to  the  Mach  cone.  The  function  f  (X)  is  shown  graphically 
in  Fig.  3.  Since  it  is  now  known  that  the  spanwise  average  of  lift  from 
tip  to  cone  is  just  .one-half  of  the  two-dimensional  value,  it  is  easy  to 
derive  formulae  for  lift  coefficient  and  center  of  pressure.  These  formu¬ 
la*  hold  also  if  the  Mach  cones  overlap  on  the  wing  without  cutting  the 
opposite  wingtip.  In  other  words,  the  formulae  hold  for  £  I  j 

where  5  is  the  mean  span,  or  average  of  the  span  at  the  leading  and 
trailing  edges. 


Cl. 


,  I  C  . 

=  l--g-  tan p.  + 


tan  S|  +  ton  Sgl 

2  J 

Sl~St 

4s 


(«) 

(8>) 


and  .Sip  are  the  spans  at  the  leading  and  trailing  ^dges  respect 
ively.  The  center  of  pressure  is  located  at  the  following  distances 
behind  the  leading  edge: 


.  .  tan»i+tanSi 

C  i-7'3  Vs  ton  P-  +  /;  Vt  2 

2"  |  -  '/2%  tan  p.  +  '/j  c/s  tan  5m- ton  St 


(9) 
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It  may  he  well  to  repeat  at  this  point  that  th6se  results  hold  only  for 
positive  rake  angles  (Si^Sfp;. 

Bent  Loading  Edge 

A  problem  of  considerable  importance  is  illustrated  in  Fig.  2. 

The  angles  and  6*,  are  not  necessarily  acute  angles,  as  in  the  case 
drawn,  hut  each  is  assumed  to  lie  in  the  range ju.  <  p  <vr-jU  .  To 
simplify  the  immediate  discussion  it  is  assumed  that  the  angle  of  the 
leading  edge  points  upstream,  so  that  6,  +  5*,  <  tT  ;  it  is  shown  later 
tnat  the  formulae  obtained  are  valid  without  this  restriction.  Evident¬ 
ly  the  wing  separates  the  problem  into  two  parts  which  may  be  treated 
independently.  Attention  may  be  confined  to  the  upper  half,  since  the 
solution  for  the  lower  half  differs  only  in  sign.  The  points  on  the 
circle  for  which  p-  (3,  ,  and  <psir-ft z  mark  tne  tangency  of  the  plane 

Mach  waves  from  the  leading  edge  with  the  Mach  cone.  By  elementary 
geometry, 

cos  /?,  =  tan  juVtan  S| 
cos  ^  =  tan  ju./'tan  £2 

From  a  consideration  of  the  "domains  of  influence",  see  (Fig.  2)  it  is 
seen  that  the  flow  in  the  regions  1^,  Ig,  III  is  as  indicated  in  the 
legend  of  Fig.  2.  As  before,  the  normal  derivative  of  w  vanishes  on 
the  wing.  All  three  components  pf  perturbation  velocity  vanish  on  the 
arc  p<-JT~  fiz  ;  on  the  other  two  arcs  the  boundary  coricLitions  are  to 

be  obtained  from  the  essentially  two-dimensional  problem  of  an  infinite 
span  airfoil  at  an  angle  of  yaw,  to  which  we  now  turn. 

Let  £  be  the  angle  between  the  main  stream  direction  and  the  lead¬ 
ing  edge  of  an  airfoil  of  infinite  span,  so  that  ~rr/z  -  $  is  the  angle 
of  yaw.  It  is  assumed  that ><«*.<  S<-rT-^c*  .  The  uniform  flow  W  may 

be  considered  as  a  superposition  of  a  uniform  flow  at  velocity  W  cos  5 
parallel  to  the  leading  ed^j,  giving  rise  to  no  perturbation,  and  a  flow 
perpendicular  to  the  leading  edge  at  velocity  W,  srW'5in3>  and  effect¬ 
ive  Mach  number  tl,  =  AJswS  •  The  effective  angle  of  attack  ) 

for  this  second  flow  is  measured  in  a  plane  perpendicular-  to  the  lead¬ 
ing  edge;  it  is  related  to  the  streamwise  angle  of  attack  (  Cf  ),  meas¬ 
ured  in  a  plane  containing  the  stream  direction  and  perpendicular  to 
the  plane  of  the  wing  for  zero  angle  of  attack,  by  the  formula 

tan  a  =  tan  cc*  sin  6 

Within  the  limits  of  validity  of  the  linear  theory  we  need  not  distinguish 
between  the  angle  of  attack  and  its  tangent,  so  that 

oc  =  ce8  sin  & 

For  a  plane  airfoil  of  infinite  span  and  not  yawed,  the  streamwise  com¬ 
ponent  (w)  of  perturbation  velocity  is 


\ij  =  ccWtcn  )i  = 


ocW 

l/  *2  r 
*  nfi  I 
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To  obtain  the  chordwlae  component  of  perturbation  velocity  for  an  air¬ 
foil  at  an  angle  of  yaw  ^/z  -  S  ,  replace  cf,Wan&  Mby  <£»,  W,  and  M, . 

The  streamwlse  component  of  perturbation  velocity  follows  from  multi¬ 


plication  by  sin  8 : 


W  = 


«.w. 


VSh 


sin 


s 


ocWsin  8 

Vm^S-I  <10J 


or 


VJ, 


1*7  = 


oo 


sm 


where,  as  before,  cos  /$  =  fan  JX/ton  $ 


The  boundary  conditions  for  the  remaining  two  arcs  of  the  upper 
semicircle  are 

w-ife-K.  ^  o«p<£, 

and 

ur  =  -^*5  =KZ  for  TT-/5,2<9<Tf 
sinP2 

The  potential  problem  which  is  now  uniquely  determined  in  the  upper  semi¬ 
circle  may  be  written  out  immediately  as  a  Fourier  series.  This  is  a 
cosine  series  only,  because  of  the  condition  that  bur/dn*0  for  <p*0  and  for 
for  <p  =YT. 


TJU  = 


K|  ^2  Kg 
TT 

/gjKi+v^  Kg 
IT 


fix 


ir 


cos  n^-^Kj /cos  nu  du  +  Kg/< cos  nu  du^ 

o  K-'fli  ^ 

+.  &/tan<r M9+0,)  _  tnn'  rsin^-/?,) 

Tf  l  l-r  cos(ip+/S,)  l-r  r,os(9-4;.) 

-  & )  ( u) 
Tt  l  i+r  cos(9+y5e)  j+r  cos(f-/S2)J 


In  Eq.  (11)  each  inverse  tangent  is  restricted  to  its  principal  values 
(~T Vt,  to  +  'trlz).  The  details  of  this  summation  are  not  much  different 
than  for  the  case  of  the  rectangular  wing,  and  need  not  be  dwelt  upon 
here. 


In  the  symmetrical  case  §|  =  the  expression  (11)  for  w  on  the 
wing  (<p»OjtT  )  simplifies  to 
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UJ 


||  ^  !  »>2  OR  1 


_2_ 

Tf 


WZL- 

sin/S 


tan 


r\ 


i-r2  cos  20 } 
fan/# 

Vi  -&/W 


(13) 


Returning  to  the  general  case  (  and  S^not  necessarily  equal)  we  seek 
the  average  value  of  w  along  the  segment  of  a  spanwise  line  (i,e.,  per¬ 
pendicular  to  the  stream)  cut  off  hy  the  Mach  cone.  It  is  necessary  to 
evaluate  integrals  such  as 


r  sin  Y 
I  -  r  cos  if 


dR 


Since  RfA-ZA/i+yi*  »  integration  "by  parts  leads  to 


I  (  -i  r  sin  V*  jrN  Tf  ,  _  .  .  T 

AJ,an  dR«-4-tanY-+(i-secr)- 


Tf-r 


where  o<V=<2 tr.  Using  this  result  it  is  found  that  the  average  value 
8 ought  is 


uj  =  -^-(l+tan  ^-sec^,)  +^(l+tan/52-sec^2)  (ls> 


Case  of  Intersecting  Envelopes 

For  treating  forward  sweep  or* dihedral,  it  is  necessary  to  discuss 
the  plane  waves  from  the  leading  edge  in  more  detail.  On  one  aide  there 
ia  a  weak  shock  wave  and  on  the  other  side  a  weak  expansion  wave;  how¬ 
ever,  in  the  linear  theory  the  distinction  "between  shock  waves  and  rare¬ 
faction  waves  disappears.*  Both  are  regarded  merely  as  surfaces  of  dis¬ 
continuity,  which  can  occur  only  across  the  envelope  of  Mach  cones  with 


*  It  is  this  fact  which  makes  possible  the  usual  two-dimensional  linear 
and  second  order  calculations,  in  which  the  pressure  is  determined  by 
local  conditions  and  does  not  depend  on  the  history  of  the  flow  up  to 
that  point.  The  entropy  increase  is  of  the  third  order  in  the  pertur¬ 
bation  velocity  components,  whereas  the  linear  theory  retains  only  the 
first  order. 
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vertex  on  the  leading  edge.  As 
the  simplest  example,  Pig.  4 
shows  a  section  perpendicular 
to  the  main  stream  for  the  case 
of  a  wing  with  dihedral  and  with 
the  leading  edge  perpendicular 
to  the  main  stream.  The  situa¬ 
tion  for  a  plane  wing  with  for¬ 
ward  sweep  would  differ  only  in 
that  the  trace  LtfiT  of  the  wing 
in  Pig.  4  would  he  straight,  and 
overlapping  of  the  plane  waves 
would  occur  over  the  bottom  arc 
as  well  as  the  top. 

The  envelope  of  all  Mach 
cones  with  vertex  on  LM  consists 
of  two  half  planes  with  traces  AC 
and  ac.  Similarly  the  Mach  cones 
with  vertex  on  MR  give  rise  to 
the  envelope  represented  by  FG 
and  fg.  It  is  important  to 
notice  that  no  arc  of  the  circle 
is  a  part  of  either  envelope. 

Within  the  bounds  of  the  linear 
theory,  shock  waves  or  rarefaction  waves  intersect  without  mutual  inter¬ 
ference,  and  the  perturbations  caused  by  each  are  additive.  In  the  region 
GBCDN  the  flow  is  uniform;  in  the  region  ABFEL  there  is  another  uniform 
flow;  in  the  region  between  FBC  and  the  circle  the  flow  is  also  uniform, 
since  the  components  of  the  perturbation  velocity  are  obtained  by  addition 
of  the  components  of  perturbation  velocity  for  the  other  two  uniform 
flows.  This  completes  the  specification  of  the  boundary  conditions  for 
the  upper  part  of  the  circle.  When  dealing  with  u,  v,  or  w,  the  function 
sought  assumes  a  constant  value  of  EF,  another  constant  value  on  CD  (bpth 
of  these  constants  obtainable  from  two-dimensional  theory},  and  a  con¬ 
stant  value  on  the  arc  PC,  namely  the  sum  of  the  other  two  constants. 

The  boundary  conditions  for  the  lower  part  of  Fig.  4  present  nothing  new; 
u,  v,  and  w  take  on  calculable  (constant)  values  on  the  arcs  Ec  and  fD, 
and  the  value  zero  on  the  arc  cf. 

Review  of  the  problem  illustrated  by  Fig.  2  now  shows  that  the  analysis 
given  holds  also  for  the  case  of  a  wing  with  forward  sweep,  that  is,  a 
wing  with  the  angle  pointing  downstream.  The  only  difference  is  a  slight 
modification  of  Fig.  2. 


In  Cott  drown  Lending  Edge 
is  Perpendicular  to  Streom 


FiG.  4 

WING  WITH  DIHEDRAL, 
SECTION  BY  PLANE 
PERPENDICULAR  TO  STREAM 
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Trapezoidal  Wing 


The  lift  coefficient  and  center  of  pressure  for 
trapezoidal  wing  shown  in  Fig„  5  may  now  V  studied, 
trailing  edges  are  perpendicular 
to  the  main  stream,  and  the  tip 
angle  ( 5  )  is  greater  than  the 
Mach  angle.  Since  the  leading 
edge  is  perpendicular  to  the 
stream,  5,  =  = ir/z 

and  the  subscript  may  he 
dropped  from  Si  and  (3t 
In  the  region  I, 


the  symmetrical 
The  leading  and 


UJ 


=  i. 


FIG.  5 


In  the  region  II,  the  average 
w  is  found  from  Eq.  (13): 


W 


l+tan/?"  secA 
2  sin/tf 


In  the  region  III, 


—  =  CSC/& 

Woo 


On  talcing  the  average  of  these  quantities,  weighted  according  to  the 
area  in  which  each  applies,  it  is  found  that 

CL  _  !  (14) 


Similarly  the  center  of  pressure  is  found  to  lie  behind  the  leading  edge 
by  the  distance 

2  = -§- (l  +  j ton  S]  (is) 

Thus  in  this  case  the  lift  coefficient  and  center  of  pressure  are  the 
same  as  if  the  wing  were  sub.lect  to  the  uniform  lift  distribution  of 
an  infinite  span  airfoil,>  The  actual  lift  is  not  uniform;  in  the  region 
I  the  lift  is  that  of  an  infinite  span  airfoil;  in  the  region  II  the 
lift  is  less,  and  in  the  region  III  the  lift  is  greater  by  just  enough 
to  compensate  for  the  decreased  lift  in  II.  As  an  example,  Fig.  6  shows 
the  apanwise  pressure  distribution  for  the  case  $=*  45*  /*,-»  3o°  s>). 
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-I  0  I 

x4  TAN  p. 


FIG.  6 

PRESSURE  DISTRIBUTION 
ALONG  A  SPAN  LINE  FOR 
TRAPEZOID  WITH  S -45°, 
M  *30°,  M  =  2 

(see  ns.  5) 


General  Symmetrical  Quadrilateral 

The  case  of  a  quadrilateral  which  is  symmetrical  about  a  diagonal, 
that  diagonal  being  parallel  to  the  stream  is  considered  here.  The  semi- 
vertex  angles  at  the  nose  and  tail,  say  §  and  respectively,  are  not 
necessarily  acute  angles  (see  Fig.  8  for  the  various  possibilities)1,  it 
is  assumed  only  that  each  lies  in  the  range  JA  <  S .<  TT*>W  . 

It  is,  of  course,  necessary  that  5*5,  <-rr . 

The  forward  pointing  triangle  is  a  special  case,  fi,  a  Tf/fc , 

It  is  also  a  special  case  of  the  trapezoid;  setting  StmQ  or  S  ~  ctnti  6 
in  Eqs.  (14)  and  (15)  leads  to 

Cl/Clk=  I 
H  =  2c/3 

Here  c  is  the  distance  from  the  vertex  to  the  trailing  edge-.  Fig,  7 
shows  the  pressure  distribution  spanwise  for  the  case'  5 
yU.  a  30°  (h  =  2)  .  Fig.  ?  Is  applicable  to  any  other  case  for  which 
tan  5=  7^  tan. yu  by  a  uniform  change  of  scale. 

For  the  more  general  quadrilateral  of  Fig,  8,  the  lift  distribution 
is  shown  by  Eq.  (12)  in  connection  with  Eq,  (3);  the  calculation  of  lift  and 


(16) 
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Aerodynamic  Characteristics  of  WingB  at  Supersonic  Speed* 

center  of  preseure  involves  a  tedious  integration.  The  most  convenient 
procedure  for  finding  the  average  pressure  over  the  wing  is  to  integrate 
along  any  line  (AB  in  Pig.  8)  parallel  to  the  trailing  edge  on  one  side. 
This  give 8  for  the  average  value  of  w,  hy  Eq.  (12) 


tf 


— 
U  = 


co 


(Z  A  -i  tan/?  JA  (\ 

JtF  tan  vpp-dtV  dt 


sm/j 


D 

where  X=t 
A 


7v 


t  = 


XsecS, 


A+  $ec/?i 

X=i  sin(juL+S,)  ^x=«ec/? 


\=^,  =  7, 


smjj. 


sin  8 
sin  (6 +6,} 


The  integration  lead*  eventually  to 

CL  2  #sin2/g-/ysin  Z$\ 

Gl00“1uj<»""7T  sin/?,  sin2^-sin/?sin2/?, 


The  angles ^3  and  t  are,  of  course,  to  he  measured  in  radians.  This 
expression  for  Ct/Ci.w  is  Symmetric  in  the  two  angles  and.  a  , 
and  therefore  is  unchanged  hy  interchanging  S  and 


The  center  of  pressure  is  at 


(l“t  cos  S|) 


(18) 


where  t  is  the  weighted  average  of  t  along  AB,  with  weights  proportional 
to  the  pressure*  It  is  found  that 


^2 

dt  +J  t  dt 


|  cos2/9>+cosfe  sin2/?  singfl-gflcosZA 
2  cos2/?,  -cos2/?  +  2sin^A  '2/?,sin20-2/tf  sin2/?, 


(19) 
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These  formulae  contain  as  special  cases  the  forward  pointing  tri- 
angle  (fi^'tr/% )  and  the  "backward  pointing  triangle  (/3=~rr/z )« 

For  these  triangles,  the  invariance  of  lift  coefficient  (and  in  this 
particular  case,  the  center  of  pressure  also)  with  respect  to  reversal 
of  the  flow  direction  may  "be  easily  verified.  It  has  already  been  shown 
that  lift  coefficient  and  center  of  pressure  for  the  forward  pointing 
triangle  are  the  same  as  if  the  pressure  were  uniform  (which  it  is  not; 
cf.  Fig.  7).  For  the  backward  pointing  triangle  the  pressure  is_  uni¬ 
form.  Reversal  of  direction  of  flow  causes  a  radical  change  in  pressure 
distribution  even  though  it  does  not  alter  the  lift  coefficient  for  these 
quadrilateral  wings. 

Another  special  case  of  interest  is  the  diamond  (^8,^  8  ;  fit  -fi)» 

In  this  case  Eqs.  (18)  and  (19)  reduce  to 


CL  sin  20- 2/^0052/^  (so) 

Cu  ITsinV? 

1  2/S  smg2/9 

Z  .  3  sin2/g-2/gcos2/g  (21) 

C  !  ~  cos  2.P 


Using  Eq§.  (20)  and  (21),  the  lift  coefficient  and  center  of  pressure 
for  a  diamond  have  been  evaluated  and  are  presented  in  the  following 
table.  It  is  seen  at  once  that  th$  property  of  invariance  under  re¬ 
versal  of  flow  direction,  which  was  found  to  hold  for  does  not  in 
general  apply  to  the  center  of  pressure.  If  such  were  the  case  the 
center  of  pressure  pf  the  diamond  would  necessarily  be  at  Z  »  c/a. 
From  Table  I  it  is  seen  that  the  center  of  pressure  of  a  diamond  act¬ 
ually  lies  forward  of  the  midpoint,  though  never  forward  of  7cJiS 
as  long  as  § 


Table  I. 


0° 

100 

200 

300 

400 

Ct/°lco 
*/•  / 

t  any  tan ^ 

.8488 

.4667 

1.000 

.8592 

.4685 

1.064 

.8720 

.4709 

1.155 

.8897 

.4743 

1.305 

50° 

60° 

70° 

80° 

90® 

°l/Cl 

*/«. 

tan Y tan// 

.9120 

.4788 

1.556 

.9376 

.4842 

2.000 

.9646 

.4905 

2.924 

.9885 

.4966 

5.759 

1.000 

.5000 

CO 

so 
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Wing  with  Dihedral 

The  case  of  a  wing  with  dihedral,  which  may  "be  combined  with  for¬ 
ward  or  backward  sweep  at  the  dihedral  point,  is  reducible  by  a  conform¬ 


al  transformation  to  the  case  of  a 
flat  wing.  In  Fig.  9,  the  angle  <p 
no  loss  of  generality).  Let  Y 
be  the  radian  measure  of  the  arc 
of  the  circle  subtended  by  the 
wing.  The  part  of  the  circle 
not  drawn  refers  to  an  independ¬ 
ent  problem  of  the  sane  type, 
with  a  different  Y  •  The 
angles  /3#  and  have  the  same 
meaning  as  in  the  plane  pase. 

In  the  case  shown,  <  / 

(which  incidentally  corresponds 
to  considerable  sweepback  i  I.’,  t hi  S 
case),  but  this  condition  is  not 
essential,  and  is  assumed  here 
merely  to  simplify  the  drawing. 

The  potential  problem  to  be 
solved  in  the/ij<p  plans  in¬ 
volves  the  following  boundary 
conditions: 


simple  bend  in  the  leading  edge  of  & 

=  O  has  been  taken  in  the  wing  (with 


FIG.  9 


y~  =  O  for  o=0  or  f=r,  0<r<l 

tu-  K,  for  O<f<-0,  ,  r «  I 

w  =  Ka  for  r=| 

w  =  0  for  r  « | 
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/  ^  7T/f° 

The  transformation  &  -  &  maps  the  relevant  part  of  the 

unit  circle  in  the  £  =  Ag*?  plane  into  the  upper  half  of  the  unit 
circle  in  the  %! -A'c*?'  plane.  It  is  clear  that  w  is  given  "by  Eq. 
(11),  with  r  replaced  "by  /£#/*'  ,  and  every  angle  (73.  t  /3j>  and  <D ) 

multiplied  "by  7T/y 

It  remains  to  specify  KtJ  in*  terms  of  the  geometry 

of  the  wing.  As  in  the  plane  case, 


COS/^i 


tan  jt 
tan  6, 


COS^g  = 


tanjx 
tan  $2 


where  S,  and  are  the  angles  from  the  main  stream  direction  to  the 
leading  edge.  The  boundary  values  and  are  given  by  the  same  ex¬ 
pression  £S  for  a  plane  wing,  except  that  the  angle  of  attack  is  now„ 
in  general,  different  for  each  plane  of  the  wing.  Writing  cL,  and 
CCZ  for  these  "local"  angles  of  attack,  one  find  that 


_  giwtgn^  „  _  ggW  tan  )X 

sin#  2"  sin/?2 

It  is  interesting  to  note  that  in  the  case  of  symmetry 

(  j  ) 

there  is  a  certain  dihedral,  namely  V*=  ZjS  ,  for  which  the  boundary 
condition  is  w  *  k  over  the  entire  arc y  ,  so  that  the  flow  iB  uniform 
in  the  whole  sector. 

The  restriction  /3,+/3 V  is  seen  to  be  non-essential  as  in  the 
previous  case  of  a  bent  leading  edge.  Also  it  should  be  pointed  out  that 
for  a  wing  with  upswept  dihedral  for  the  upper  surface)  the, 

lift  is  decreased  in  magnitude. 

An  application  of  these  formulae  may  be  found  in  the  perpendicular 
vane  at  the  tip  of  a  rectangular  wing,  the  vane  being  large  enough  to 
project  through  the  Mach  cone.  This  may  be  regarded  as  an  example  of 
a  wing  with  dihedral,  one  of  the  angles  of  attack  being  zero.  If  the 
vane  extends  both  above  and  below  the  wing,,  the  lift  remains  constant 
out  to  the  end  of  the  wing.  If  the  vane  is  confined  to  either  the  top 
or  the  bottom  of  the  wing,  the  lift  decreases  moving  along  a  span  line 
toward  the  wingtip;  the  lift  at  the  tip  is  1/3  of  the  lift  in  the  cen¬ 
tral  region  of  the  wing,  and  the  spanwiBe  average  from  tip  to  Mach  cone 
is  found  by  integration  to  be  /  /  %\  of  the  lift  in  the  central 

region.  '  373/ 
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I.  ABE0M5AMIC  CHARACTERISTICS  01  SOLID  RECTANGULAR  AIREOILS  AT  SUPERSONIC  SPEEDS* 

by  E.  A.  Bouncy 


The  lift,  drag,  moment,  and  center?of-presBure  characteri*tic«  of 
rectangular  airfoils  can  readily  he  obtained  for  angles  of  attack  below 
that  at  which  the  shock  wave  detaches  from  the  leading  edge  of  the  air¬ 
foil,  This  report  presents  expressions  for  the  above  characteristics 
and  establishes  optimum  values  of  aspect  ratio,  angle  of  attack,  and 
lift-drag  ratio  for  any  given  conditions  of  allowable  stress,  airfoil 
cross-sectional  shape  and  Mach  number  as  well  as  method  of  supporting 
the  wing,  i.e.  by  the  entire  base  or  on  a  hub. 

As sumptions 

It  was  necessary  to  make  certain  simplifying  assumptions  to  keep 
the  various  axpretsions  reasonably  simple  and  yet  accurate.  They  are 
as  follows: 


1.  The  change  in  Mach  angle  with  positive  angle  of  attack 
over  an  airfoil  of  finite  thickness  and  aspect  ratio  was  ignored.  The 
angle  will  be  higher  on  the  lower  surface  and  lower  on  the  upper  sur¬ 
face,  thereby  offsetting  each  other  to  a  great  extent. 

2.  The  possibility  of  secondary  tip  effects  originating  at 
the  point  of  maximum  thickness  of  a  double  wedge  airfoil  for  example 
was  not  considered. 

3.  .Consideration  of  the  phenomena  of  separation  necessarily 
was  omitted  due  to  the  lack  of  knowledge  of  this  effect  in  supersonic 
flow.  This  factor  can  cause  the  center  of  pressure  expressions  ob¬ 
tained  herein  to  be  somewhat  in  error,  particularly  at  high  angles  of 
attack.  A  study  of  the  effects  of  separation  is  contained  in  reference 
(12). 


4.  A  constant  skin  friction  drag  coefficient  of  =*  .00265 
per  square  foot  of  wetted  area  is  used  throughout.  This  is  the  value 
obtained  by  von  Earman  for  a  Reynolds  number  of  about  20,000,000.  It 
is  recognized  that  the  actual  value  may  be  considerably  different  than 
this  and,  inasmuch  as  skin  friction  is  the  predominant  factor  in  the 
drag  of  thin  airfoils,  can  cause  a  corresponding  difference  in  the 
values  for  optimum  conditions. 

*  This  paper  is  a  revision  of  CM-247  which  was  originally  published  as 
an  internal  memorandum  of  the  Applied  Physics  Laboratory.  Slight  re¬ 
vision  was  necessary  to  complete  the  discussion  of  the  subject. 
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5.  The  expression  for  aspect  ratio  correction  to  lift  curve 
slope  for  an  infinitely  thin  flat  plate  was  used  in  deriving  optimum 
conditions  for  simplification.  The  maximum  error  that  could  result  from 
this  assumption  is  about  3  per  cent  for  practical  airfoil  sizes, 

6.  Throughout  the  analysis  the  assumption  is  made  that 

Sin  CC  =  GC  and 
Cos  0C  =  I 

This  will  cause  no  sizable  error  for  angles  of  attack  up 

to  15  degrees. 


Nomenclature 

The  following  nomenclature  is  used  throughout  this  paper: 
b  =  wing  semi- span 
c  =  chord 

ft  =  2b  *  aspect  ratio 
c 

M  =  Mach  number 
44-=  Vm2  -  1 
R.  =  ft  x  44- 

M  -  Mach  angle  =  Sin-1  I  =  Tan-1  i 

M  -M- 

p  =  local  static  pressure  at  any  point  on  the  airfoil 
P0  =  free  stream  static  pressure 

=  free  stream  density 

V  =  velocity  in  ft/sec. 

<1  “  oVS  =  Y  p0  M§  *  free  stream  dynamic  pressure 

2  “2 
L  »  lift 

A  =  2bc  =  wing  area 

-  L  =  lift  coefficient  of  finite  span  airfoil 
qA 


Aerodynamic  Characteristics  of  Solid. 

Rectangular  Airfoils  at  Supersonic  Speeds 

*  lift  coefficient  of  an  infinite  span  airfoil 
D  =  drag 

C5  =  _D  r  drag  coefficient  of  finite  span  airfoil 

qA 

M  =  moment  about  leading  edge  of  wing 

C#  =  M  =  moment  coefficient  about  leading  edge 
cqA 

c.p.  =  center  of  pressure  measured  from  leading  edge 
Gp^  *  skin  friction  drag  coefficient 
Cp*  =  form  drag  coefficient 
Cp^  =  wave  drag  coefficient  due  to  lift 
X  ^  =  coefficient  for  taper  in  thickness 
K-^  -  coefficient 'for  wing  form  drag 
K4  *  coefficient  for  wing  strength 
w  =  Op  q  =  wing  loading 
K  3KxKiX4v 

- S5T 

S  =  stress 

s  =  section  modulus 

ti  *  wing  thickness  at  root 

tg  -  wing  thickness  at  tip 

&  -  angle  of  attack 

ft  3  semi-vertex  angle  of  leading  edge 

0  =  local  angle  between  any  point  on  the  surface  of  the 
airfoil  and  the  free  stream  direction 

Y  =  ratio  of  specific  heats 

p/=  P~Po  3  AP.  =  pressure  coefficient 

”  <T  q 

Subscripts 

U  =  upper  surface  of  airfoil 
L  =  Lower  surface 
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F  -  forward  half  of  douole  wedge  airfoil 

R  »  rear  half 
o  =  free  stream  conditions 


Discussion 


A.  Basic  equations  of  lift,  drag,  moment,  and  center  of  x>res~ 
sure  for  infinitely  thin  flat  plates  and  for  airfoils,  both  with  infinite 
and  finite  aspect  ratios. 


Various  methods  (references  4,  10,  and  11)  have  been  developed  for 
determining  the  pressure  on  an  airfoil  in  supersonic  flow  where  a  shock 
wave  has  formed  and  is  attached  to  the  leading  edge.  However  they  do 
not  all  lend  themselves  to  convenient  handling  for  purposes  of  develop¬ 
ing  expressions  for  lift  and  drag  of  the  airfoils.  The  one  exception 
to  this  is  the  Busemann  "second  order  approximation"  (see  ref.  4  and 
10)  which  is  usually  written  as  follows: 

fypo=  I  +  23f  cosec  2>it>8+^sec2jx0  cosec2  2jlio(#+cos22>i)02  (D 


Transforming  "by  use  of  the  fundamental  relationships, 

sin  >o=^  tan  6  =  !^  C08®=^M0  and 

a  convenient  expression  for  &P  is  obtained  as  follows: 

<1 


Ap  _  2  t  VM4+(M2-2)g  02 

q.  vpzj  2(M8-0* 


®  C,9  +C£©2 


(2) 


This  method  is  found  to  be  very  accurate  when  compared  to  thv»  theo¬ 
retically  correct  method  of  "patching  curves"  of  reference  (4)  two  dimen¬ 
sional  method  of  characteristics)  with  the  error  approaching  about  -2 
per  cent  in  lift,  drag,  and  moment  at  about  60  per  cent  of  the  detachment 
angle.  Above  this  the  error  increases  roughly  in  a  parabolic  manner  un¬ 
til  at  the  detachment  angle  it  may  amount  to  -10  to  -13  per  cent,  depend¬ 
ing  on  the  geometry  of  the  airfoil.  The  angle  of  attack  at  which  detach¬ 
ment  occurs  (called  detachment  angle)  decreases  with  decreasing  Mach 
number  and  increasing  leading  edge  wedge  angle  and  may  be  determined 
for  any  condition  from  reference  (4).  The  error  in  center  of  pressure 
location,  neglecting  separation,  when  compared  to  the  accurate  method 
is  negligible  at  any  angle  of  attack  up  to  detachment. 

From  this  expression,  it  is  possible  to  determine  the  aerodynamic 
characteristics  of  any  shape  of  airfoil  of  infinite  aspect  ratio 
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and  if  the  loss  in  pressure  at  the  tips  for  finite  aspect  ratio  wings 
is  also  known,  then  they  can  he  found  for  any  airfoil  of  any  aspect 
ratio,  For  example,  for  any  wing  of  infinite  aspect  ratio  and  symmet¬ 
rical  about  its  chord  line,  the  lift  coefficient  can  he  determined 
simplj'-  by  the  inspection  of  equation  (2)  as  being 

Cl=2C|QC=^T 

(8  a  ±<L±/3  depending  on  surface) 

The  lotfs  in  pressure  at  the  tips  due  to  flow  from  the  lower  to 
the  upper  surface  has  been  determined  in  references  (2)  and  (8)  and  is 
shown,  in  Fig.  1.  It  can  be  seen  from  this  figure  that  the  assumption 
of  a  linear  loss  in  lift  from  the  limit  of  the  Mach  angle  to  the  tip 
is  justified  and  that  the  lift  inside  the  Mach  cone  will  therefore  be 
one- half  of  the  two  dimensional  value  for  the  same  area. 

With  conditions  of  pressure  known  at  every  point  on  the  airfoil, 
it  is  now  possible  to  develop  expressions  for  lift,  drag,  moment  and 
center  of  pressure  for  any  type  of  airfoil.  These  expressions  are 
listed  in  Table  I  and  are  good  for  all  airfoils  that  are  symmetrical 
about  both  the  chord  line  and  a  line  perpendicular  to  the  chord  at  the 
midpoint.  Sample  derivation  of  some  of  these  expressions  will  be  found 
in  the  section  on  "Derivations". 

B, ..Optimum  desi/rc  conditions  for  wings  of  finite  thickness 
and  aspect  ratio  for  (a)  wings  supported  over  their  entire  base  and 
(b)  wings  supported  by  a  hub. 

In  determining  the  optimum  design  conditions  of  aspect  ratio,  drag, 
and  angle  of  attack  for  various  values  of  allowable  wing  stress  apd  type 
of  airfoil  (which  establish  the  value  of  K3),  and  Mach  number,  the  meth¬ 
ods  of  references  (7)  and  (9)  were  used  for  the  most  part. 

It  should  be  pointed  out  that  the  problem  of  establishing  optimum 
conditions  is,  in  this  paper,  principally  of  academic  interest  for  two 
reasons.  First  the  optimum  type  of  airfoil  will  not  be  solid,  but  rather 
of  monocoque  conotruction  due  to  the  excessive  weight  of  a  solid  wing 
and  secondly,  the  optimum  angle  of  attack  upon  which  all  of  the  other 
conditions  depend  is  different  for  a  wing  on  a  body  than  for  a  wing 
alone.  Various  trends  can  be  studied  however,  which  will  be  applicable 
regardless  of  these  factors. 


Derivations 


A,  In  this  section,  the  method  of  derivation  of  the  express¬ 
ions  given  in  Table  I  for  lift  coefficient  and  center  of  pressure  of  an 
airfoil  of  finite  thickness  and  aspect  ratio  will  be  shown.  Several 
of  the  other  expressions  will  be  found  in  the  various  references  (see 
references  I,  2,  3,  and  10)  although  perhaps  not  in  exactly  the  same 
form.  For  the  sake  of  brevity,  the  double  wedge  section  will  be  used 
in  the  following  illustrations  and  the  final  expressions  generalised 
by  a  method  suggested  by  Busemann. 


1 


if 
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e3*  cc  +/? 

Q*  -  0.-/8 


Therefore,  the  two  dimensional  pressure  difference  over  the  front  and 
rear  halves  of  the  airfoil, 

p®  P©=  “l  «*«*«/* 

i  (4) 

p4'p4=2Ci<r"4C2°^ 

Remembering  that  the  lift  in  region  0  is  equal  to  half  of  that  in  region 
A  (for  a  like  area)  and  likewise  the  lift  in  region  D  is  half  of  that  of 
B,  thsn 

Ul=  2bc 


(5) 
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uaing  the  expressions, 

C=H  and  R= 

th9  coefficient  of  lift  "becomes 

CL=  2C|0c 


(6) 


(7) 


This  expression  applies  to  a  double  wedge  airfoil  only.  However,  Buse- 
mann  has  pointed  out  that  the  expression  can  he  made  general  for  all  air¬ 
foils  that  are  symmetrical  about  their  chord  line  and  a  line  perpendicular 
to  the  chord  line  at  the  midpoint  by  substituting  the  parameter  of  cross- 
sectional  area  divided  by  the  chord  squared  in  place  of  the  half  wedge 
angle.  Therefore  for  a  double-wedge  airfoil, 


Acs 

C2 

/S 


tc  _  J_  _0_ 
x?  -  2c  =  2 


(8) 


Values  of  A‘  for  various  types  of  airfoils  are  noted  in  Table  I 


Therefore  if 


o  2Cg 

c.=*and  C3=-£p 


(9) 


Inasmuch  as  the  lift  coefficient  for  infinite  aspect  ratio  is 

C.  =  —  (10) 

-M- 

then  the  aspect  ratio  correction  to  2  dimensional  lift  (and  drag  due 
to  lift;  is 

CT=  '“2R(f_C*A)  (u) 

u00 


The  center  of  pressure  is  obtained  by  summating  the  center  of  pres¬ 
sure  of  each  region  times  the  total  pressure  in  that  region  for  the  en¬ 
tire  airfoil  and  dividing  by  the  total  pressure  or  lift  force. 
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Referring  again  to  the  sketch  of  the  airfoil, 

.  v  /  A  L  t  jk  !_  2,C  \ 


m  _j_/2b^4b-^\ 

'  c  4  6  \2b  +  2b“/ 

(2JL)  -i,i 

\  c  4  2  6  ^2b-^  +  2b-§ 


3R~2 

I2R-6 


9R-I4 

I2R-I8 


fc.p.\  2. 

\  c  /©"  3X  2  3 

fc-P\  J_,  5.  1  1_ 

'  c  ©  2  9  X  2  “  9 


Area  of  each  region 


A  A  jjQ 

®  =  (2b“SS:1¥  =  (2R-I)2r’ 


^2b-§>H2R-3>i 


(T)  _  JL  -  SL  - 
^  ~  2M  2  4M-  2R 


(H)  -  iL  -  3cf  ^  J*L 
^  ~  2f*  2  4M-  ”  32R 

Therefore  the  center  of  pressure  will  he,  (omitting  factor  he 
from  top  and  bottom^  2R 


^p.J2C,c+4C2|yf)|j'^)(2R-|)4  xi^t(2C|Ct -4C2a/5)[(^lA)(2R-3)+|  x  1  xs] 
(2C,oc+4C2k£)  |2R-I  + 1 X |J  +  (2C|ff -  4C zoc/3)  j2R  -3  +  ~  x  si 
=  s-3  -c*a'(r-i) 

2R  - 1  +  C3  A'  (W) 

This  expression  likewise  is  general  for  all  symmetrical  airfoils 
when  expressed  as  a  function  of  A’  rather  than^S  , 
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B.  For  a  given  value  of  allowable  stress  at  the  root  section, 
there  will  he  an  optimum  aspect  ratio  which,  for  a  given  lift,  will  pro¬ 
duce  the  least  drag.  At  low  aspect  ratios,  the  tip  loss  will  he  great 
and  hence  the  area  for  a  given  lift  will  he  large  while  at  high  aspect 
ratios  the  thickness  and  hence  the  drag  for  a  given  stress  will  he  large; 
therefore  the  optimum  aspect  ratio  will  he  somewhere  in  between  these 
two  conditions  and  is  obtained  as  follows: 


D  = 


(15) 


(<jC  in  degrees  from  this  point  on.) 
Wow  for  a  given  lift, 


L  =  Cl*  ClA®  55  CL 

Substituting  into  (15) 


(16) 


a  = 


A 


(**>.) 


(17) 


=  .00122  g2(ft-/CLn)q.Aco  Co,o 

M(VcJ  *<VO  (CL/ClJ  da) 


and  D  .00)22  a:2 

"  *  4ClAJ 

Now  for  a  double  wedge  airfoil  the  section  modulus  is 

S  a  *  m.  (wbc)b 

24  S  25 


(19) 


(20) 


where  the  spanwise  center  of  pressure  is  conservatively  assumed  to  he 
at  one-half  the  distance  to  the  tip.  From  this  expression  the  thickness 
ratio  becomes  a  » 

o  3wRe  K4wRz 


(21) 


where  is  a  function  of  the  cross-sectional  shape  and  is  given  in  Table 
I  for  various  airfoils. 

Equation  (19)  now  becomes, 

D  _  .00122  oc2  K3Rc  +CDfM  (22) 

q.A«"  m-(Cl/clJ 

NOTE 

’AjES  32  &  33 
transposed 
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The  parameter  Is  plotted  in  Figs.  12  and  18  to  show  the  variation 

in  thickness  ratio. 


From  the  expressions  for  drag  coefficient 

sARb 

w~+c°f 


r  .00122a2  |  ' 

Cd-  n  (l-gR)  + 


W- 


and  lift  coefficient 


„  .0698oc  I  , 

CL  =  — it—  (  ~  2R 


-M- 

the  drag-lift  ratio  "becomes' 


j-  =  ,01745  a  + 


KjAR^  +  CofM- 


.0698(1“  gp)a 


(28) 


Differentiating  with  respect  to  CC  and  optimizing,  the  optimum  angle 
of  attack  "becomes, 


OC0pt.  *  28.6 


2R(K3/«b+CDfM’ 


2R  -I 


(29; 


and  the  minimum  drag- lift  ratio  "becomes 

D 


(U\  _  ,  ,  2R(K5«Rb+CofM-) 

vL^min"-50V  2R-I 


K3/«"+  Cp,*f _ 

?  q(2r-1^  i/gR(K3ARb+Cp^4j 
'  V2R  '*  2R_| 


and 

(~)  = 
V  n  fUAY 


2R(K3ARD+CDfM) 

2R-I 


(30) 


2R-I 


D  'MAX  V  2R(K3ARb+C0fftH 


(31) 


The  interesting  relation  "between  maximum  lift-drag  ratio  and  optimum 
angle  of  attack  is  to  "be  noted. 

aopt  <deg}  = -—•  or  aopt  (rod.)  =  — - 

(  /d)max  (  ' 

This  relation  will  "be  true  for  a  wing  alone  or  for  a  wing  and  "body  eom- 
"bination. 


From  the  above  expressions,  the  maximum  li^t-drag  ratio  and  optimum 
angle  of  attack  are  found  to  "be  very  nearly  alike  (within  ±  3  per  cent) 
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where  K3  = 


K^Kj^w 
1  S 


(23) 


By  a  more  arduous  means,  it  may  be  shown  that  for  a  wing  supported 
by  a  hub  or  shaft,  the  exponent  of  the  aspect  ratio  will  be  approximately 
3/2  for  values  of  shaft  diameter  equal  to  1  1/2  -  2  times  the  maximum 
thickness  at  the  root.  This  is  an  approximation  (because  the  value  must 
again  be  2  for  hub  diameters  equal  to  or  greater  than  the  chord  length), 
but  it  will  be  shown  that  the  method  of  support  does  not  affect  the  aero¬ 
dynamic  values  to  any  great  extent,  but  simply  changes  the  optimum  aspect 
and  tnickness  ratios. 


Optimizing  equation  (22)  for  aspect  ratio, 

d  /Kj^  +  CofMs 

d*W(C'-/CLJ  )  =  ° 

Using  the  expression  for  aspect  ratio  correction  to  lift  curve  slope 
of  a  flat  plate  for  values  of  R>1,  the  expression  for  optimum  aspect 
ratio  for  a  wing  supported  over  the  entire  base,  becomes, 

W-3=^  RZ(4R-3)  (24) 

and  fos  wings  supported  by  a  hub , 


R^(3R-f-) 


(25) 


These  expressions  are  plotted  in  Figs.  7  and  13  to  show  the  varia¬ 
tion  of  optimum  aspect  ratio  with  Mach  number  and  the  aerodynamic-strength 
parameter  K^.  The  corresponding  expressions  for  values  of  R  <  1  will  not 
be  derived  here  because,  for  the  most  part,  such  values  are  outside  of  the 
realm  of  practical  aspect  ratio  -  Mach  number  combinations. 


From  the  development  of  optimum  aspect  ratio,  it  becomes  apparent 
that,  for  wings  supported  over  their  entire  base,  that 

K,r2=  k3/«2  ’  <*> 


and  for  wings  supported  by  a  hub, 

K,72  =  K3/R%* 


(27) 
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for  "both  methods  of  support  and  the  mean  values  are  Shown  plotted  in 
Figs.  9,  10,  15,  and  16. 

In  order  to  determine  the  value  of  wing  loading,  w  to  be  U9ed  in 
evaluating  Kg,  it  Is  necessary  to  know  either  the  weight  and  wing  area 
or  the  lift  coefficient  and  dynamic  pressure.  The  dynamic  pressure  will 
depend  on  the.  choice  of  altitude  and  Mach  number.  Inasmuch  as  the  opti¬ 
mum  aspect  ratio  and  angle  of  attack  have  been  established,  the  optimum 
lift  coefficient  can  readily  be  found  and  develops  to  be  very  nearly  in¬ 
dependent  of  .the  method  of  support  and  the  parameter  Kg.  The  mean  value 
is  plotted  in  Figs.  8  and  14. 


Results 

Expressions  for  lift,  drag  and  moment  coefficients  and  center  of 
pressure  for  thin  flat  plates  and  airfoils  of  infinite  and  finite  aspect 
ratio  are  given  in  table  I.  The  lift  curve  slope  correction  and  center 
of  pressure  expressions  for  a  flat  plate  will  be  found  plotted  in  Figs. 

3  and  4.  The  effect  of  thickness  on  lift  has  not  been  added  due  to  its 
small  magnitude  and  the  fact  that  it  does  not  lend  itself  to  a  general 
expression  for  plotting,  however  the  effect  on  center  of  pressure  is  quite 
sizable  and  two  sample  cases  have  been  plotted  in  Fig.  4.  The  express¬ 
ions  of  table  I  are  exact  (within  limits  of  accuracy  of  the  Busemann  sec¬ 
ond  order  approximation)  down  to  a  value  of  R  -  1.  Below  this  value,  a 
linear  extension  to  0  will  be  approximately  correct. 

The  loss  in  lift  at  the  tips  due  to  flow  around  the  tips  as  found 
in  references  (2)  and  (8)  and  is  plotted  in  Fig.  1.  A  straight-line 
variation  is  also  shown  to  illustrate  the  fact  that  the  lift  in  the  af¬ 
fected  region  (insJ.de  the  Mach  angle  from  the  tip)  is  one-half  of  the 
amount  for  two  dimensional  flow  for  the  same  area. 

'A  correction  factor,  K\  for  tapered  (in  thickness)  airfoils  of 
rectangular  planform  to  be  used  in  the  expression  for  Kg  is  given  in 
Fig.  5  Sad  the  terms  explained  in  Fig.  6.  These  two  curves  were  taken 
directly  from  reference  (6). 

Figures  7  to  12  and  13  to  18  represent  the  optimum  aerodynamic  char¬ 
acteristics  of  symmetrical  wings  as  a  function  of  Mach  number  and  Kgl( 
the  aerodynamic-strength  parameter,  for  wings  supported  over  their  en¬ 
tire  base  and  for  wings  supported  by  a  hub.  The  similarity  in  magnitude 
of  the  aerodynamic  coefficients  for  wings  supported  over  their  entire 
base  as  compared  to  wings  supported  by  a  hub  is  to  be  noted.  Optimum 
aspect  ratio  and  thickness  ratio  are  different  for  the  two  cases,  but 
compensate  for  each  other  so  that  the  coefficients  are  nearly  the  same. 

Conclusions 


The  Busemann  second  order  approximation  theory  for  the  pressure 
over  a  two  dimensional  wing  and  the  known  theory  for  pressure  loss  near 
the  tips  of  a  finite  span  airfoil  provide  a  convenient  method  for  finding 
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the  aerodynamic  coefficients  for  flat  plates  and  airfoils  of  finite  and 
infinite  aspect  ratio  at  angLes  of  attack  below  the  detachment  angle  in 
supersonic  flow.  This  second  order  theory  is  very  accurate  for  angles 
of  attack  up  to  about  sixty  per  cent  of  the  detachment  angle  after  which 
it  departs  from  the  exact  theory  at  a  more-or-less  parabolic  rate. 

By  a  method  also  suggested  by  Busamann,  it  is  possiole  to  general¬ 
ize  the  expressions  to  include  any  type  of  airfoil  section  which  is 
symmetrical  about  its  chord  line  and  a  line  perpendicular  to  the  chord 
line  at  its  midpoint. 

Decreasing  the  aspect  ratio  will  decrease  the  lift,  wave  drag  (here¬ 
in  defined  as  the  drag  due  to  lift),  and  moment  coefficients  and  cause 
the  center  of  pressure  to  move  forward.  Increasing  the  thickness  will 
increase  the  lift  and  drag  coefficients  of  airfoils  of  finite  aspect 
ratio  very  slightly,  but  will  decrease  the  moment  coefficient  and  cause 
the  center  of  pressure  to  move  forward  quite  markedly  for  airfoils  of 
finite  and  infinite  aspect  ratio.  Increasing  Mach  number  will  either 
increase  or  decrease  the  coefficients  depending  on  the  magnitude  of  the 
Mach  numbers  being  considered  and  the  aspect  ratio. 

The  expressions  for  lift  check  very  well  with  available  test  data. 
The  center  of  pressure  expressions,  however,  make  no  allowance  for  sep¬ 
aration  effects  (and  other  minor  factors  noted  in  items  1  and  2  under 
’’Assumptions1’)  and  therefore  will  be  somewhat  in  error,  the  error  being 
a  function  of  the  angle  of  attack.  The  expressions  will  make  it  possible 
to  determine  the  magnitude  of  the  separation  effects  on  center  of  pres¬ 
sure  location  however. 

Using  the  information,  developed  above,  the  optimum  conditions  of 
aspect  ratio  and  angle  of  attack  and  the  corresponding  coefficients  for 
rectangular  wings  supported  by  the  entire  base  and  supported  by  a  hub 
can  be  determined.  For  a  given  value  of  the  aerodynamic-strength  para¬ 
meter  £3,  the  aerodynamic  coefficients  are  not  a  function  of  the  method 
of  support,  being  almost  identical  for  either  method.  The  aspect  ratio, 
however,  will  be  higher  for  a  wing  supported  by  a  hub,  the  thickness 
ratio  also  being  higher  in  general  to  compensate  for  this  effect  stress- 
wise.  Other  effects  may  be  noted  by  a  study  of  Fig.  7  through  19. 

The  expressions  for  lift,  drag  and  moment  coefficients  and  center 
of  pressure  will  be  found  in  Table  I.  Note  that  an  increment  for  slrln 
friction  has  been  added  to  the  drag  equations  which  will  not  appear  in 
the  derivation  of  the  expression. 
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III.  THE  REVERSE  DELTA  PLANFORM 


A  section  on  the  "reverse  delta"  type  of  wing  as  discussed  in 
CM- ^58  was  to  have  been  included  in  this  report,  however  the  assumpt¬ 
ions  made  in  that  paper  were  not  exactly  correct,  and  therefore  only 
a  short  discussion  of  this  planform  will  he  included  here. 

The  reverse  delta  wing  has  the  planform  of  an  isoceles  triangle 
with  the  base  facing  into  the  airstream.  The  angle  of  the  rear  sur¬ 
face  is  greater  than  the  Mach  angle  so  that  no  tip  loss  effects  will 
be  felt  over  the  surface  of  the  wing.  It  was  assumed  in  the  refer¬ 
ence  report  that  thiB  therefore  constituted  a  wing  whereon  only  two- 
dimensional  flow  existed  and  the  characteristics  corresponded  to  the 
two-dimensional  case.  Actually  however,  the  pressure  on  a  surface 
behind  a  yawed  corner,  such  as  the  aft  portion  of  this  type  of  wing 
when  constructed  with  a  double-wedge  section  (region  B  of  figure)  is 
a  function  of  the  component  of  Mach  number  which  is  normal  to  the 
corner,  Mjj.  Furthermore,  the  pressure  inside  the  Mach  cone  created 
at  the  point  of  maximum  thickness  (region  C  of  figure)  is  of  the 
conical  flow  type  and  must  be  computed  by  the  conical  field  method. 
Only  in  region  A  is  the  flow  of  pure  two-dimensional  character. 

Because  of  the  various  types  of  flow  it  will  be  difficult  to 
develop  a  general  expression  for  form  drag  and  moment  coefficient. 

The  design  should  not  be  overlooked  however,  inasmuch  as  it  has  dis¬ 
tinct  structural  advantages  and  the  drag  may  be  very  little  worse  than 
the  two-dimensional  case.  The  lift  coefficient  will  be  very  close  to 
the  two-dimensional  value,  being  equal  to  it  for  zero  thickness, 

(See  figure  below) 


1 


\ 


\ 


\ 


\ 


\ 


*  "Aerodynamic  Characteristics  of  Reverse  Arrow  Vings  at  Supersonic 
Speeds",  APL/JHC  (May  1946). 
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IV.  LIFT  AND  DBAS  CHARACTERISTICS  OF  DELTA  NINOS  AT  SUPERSONIC  SPEEDS 


by  E.  A.  Bonney 


Abstract 


The  lift  characteristics  of  the  delta  type  of  wing  with  the  leading 
edge  both  inside  and  outside  of  the  Mach  cone  from  the  vertex  of  the  wing 
have  been  determined  by  Stewart  (l)»  Snow  (2),  and  others,  and  the  drag 
characteristics  (for  double  wedge  section  only)  have  been  determined  by 
Puckett  (3).  Using  the  methods  of  references  (4)  and  (5),  this  informa¬ 
tion  is  herein  used  to  determine  optimum  conditions  of  lift-drag  ratio, 
angle  of  attack,  etc.,  for  comparison  with  other  wing  planform  shapes. 

The  actual  pressure  distribution  for  wings  of  finite  thickness  ia  not 
known  at  this  writing,  and  therefore  moment  and  center  of  pressure  data 
are  not  included  in  this  report.  It  is  shown  that,  except  for  relative¬ 
ly  large  thickness  ratios  at  low  Mach  numbers,  the  delta  type  of  wing 
has  lower  values  of  maximum  lift- drag  ratio  than  the  reverse-arrow  type 
for  which  the  flow  is  entirely  two-dimensional.  The  relatively  large 
root  chord  length  required  for  a  given  lift  is  another  disadvantage  of 
the  delta  wing. 

The  principal  advantage  of  this  type  of  airfoil  is  in  its  use  as  a 
tail  surface  due  to  the  aft  center  of  pressure  location. 

In  two  reports  (6)  (7)  concerning  the  lift  of  delta  wings,  it  has 
been  mentioned  that  the  resultant  force  on  the  surface  of  delta  wings 
entirely  inside  the  Mach  cone  will  be  tilted  forward  of  the  normal  to 
the  chord  of  t.ne  airfoil  by  an  amount  which  depends  on  the  ratio  of 
the  complement  of  the  sweepback  angle  to  the  Mach  angle,  approaching 
a  limit  ofcC/2  ahead  of  the  normal  for  a  sweepback  angle  of  90  degrees. 
This  is  due  to  the  subsonic  effect  wherein  the  normal  Mach  number  con¬ 
trolling  the  pressure  is  less  than  1,  and  will  only  be  possible  where 
the  leading  edge  of  the  wing  is  rounded  to  permit  suction  peaks  simi¬ 
lar  to  the  subsonic  case. 

This  effect  has  not  been  realised  in  any  tests  to  date  however, 
and  the  following  analysis  assumes  that  the  resultant  vector  was  al¬ 
ways  at  right  angles  to  the  chord  of  the  wing. 
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Assumptions 


The  assumptions  noted  in  reference  (4)  concerning  separation,  skin 
*  friction  drag  coefficient  and  linearizing  trigonometric  functions  will 
also  he  applicable  in  this  report  and  in  addition: 

1.  The  lift  curve  slope  correction  for  delta  airfoils  is 
assumed  to  he  independent  of  thickness  ratio.  No  sizeable  error  will 
result  from  this  assumption  because  of  the  small  thickness  ratios  which 
are  used. 


2.  The  slope  of  the  lift  curve  as  derived  by  Stewart  is  at 
an  angle  of  attack  of  zero  degrees.  It  is  assumed  here  that  the  slope 
is  constant  with  angle  of  attack.  This  assumption  should  introduce  no 
sizeable  error  at  the  angles  for  maximum  lift-drag  ratio  which  are  con¬ 
sidered  herein. 


Nomenclature 

The  nomenclature  of  references  (4)  and  (5)  is  used  throughout  with 
the  following  additions  (See  also  Fig.  13): 

0  =  angle  between  the  leading  edge  and  a  normal  to  the  flow 
direction 

n  -  Tan  0  =  parameter  of  leading  edge  sweepback  angle  compared 
44-  with  Mach  angle 

£J0=  planform  semi-vertex  angle  =90-0 

c  =  basic  chordwise  dimension  measured  from  the  apex  to  the 
position  of  the  tip 

1  =  root  chord 

r  =  chordwise  distance  from  most  rearward  point  to  point  of 
maximum  thickness  divided  by  c 

a  =  l.c_-_ll 

c 

b  =  (l-r)c 
•  1 

"C  =  thickness  ratio  =  t/c 


y  =  distance  from  wing  center  line  to  center  of  pressure  of 
the  semi-span. 
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Discussion 


The  delta  airfoil  consists  of  a  planform  shape  which  is  symmetrical 
about  it*  centerline,  pointed  at  the  front,  has  a  swepthack  leading  edge 
and  straight  trailing  edge  normal  to  the  center  line. 

For  the  case  where  the  leading  edge  is  outside  of  the  Mach  cone, 
the  lift  remains  identical  in  value,  hut  not  in  distribution  ( 2 )  to  the 
two-dimensional  case.  The  form  drag,  however,  is  slightly  higher,  hence 
the  maximum  lift-drag  ratio  will  be  slightly  lower  than  for  the  reverse 
arrow  type  for  which  two-dimensional  lift  and  drag  apply. 

For  delta  shapes  where  the  leading  edge  is  inside  the  Mach  cone, 
the  lift  curve  slope  will  decrease  from  the  two-dimensional  value  by 
an  amount  which  is  a  function  of  the  ratio  of  the  leading  edge  sweep- 
back  angle  to  the  complement  of  the  Mach  angle  (defined  by  the  parameter 
nnB)  as  shown  in  Fig.  1  (1).  (Note  that  when  n  =  1,  the  Bweepback  angle 
it  just  equal  to  the  complement  of  the  Mach  angle,  and  n  =  0  corresponds 
to  a  straight  leading  edge.;  The  form  drag  will  also  be  a  function  of 
"n*  as  well  as  being  a  function  of  the  thickness  ratio,  location  of  point 
of  maximum  thickness,  and  Mach  number.  The  coefficient  of  form  drag  was 
determ  ned  in  reference  (3)  and  iB  shown  in  Fig.  2.  It  can  be  sesn  that 
for  la_ge  values  of  n,  the  form  drag  is  greatly  reduced  from  the  optimum 
two-dimensional  value,  where  n  =  0.  Hence,  for  conditions  of  high  wing 
leadings  where  the  required  thickness  ratio  becomes  large  and  form  drag 
Oecomes  a  big  proportion  Of  the  total  drag  for  optimum  lift-drag  condi¬ 
tions,  the  delta  type  of  wing  design  will  theoretically  Bhow  up  to  ad¬ 
vantage  over  the  reverse  arrow  type  in  the  lower  range  of  Mach  numbers, 
however,  for  most  practical  thickness  ratios  and  Mach  numbers,  the  re¬ 
verse  arrow  will  still  give  higher  lift-drag  ratios  as  shown  in  Figs.  3, 

4,  and  7.  A  definite  disadvantage  of  the  delta  wing  where  n  is  large  is 
the  relatively  large  root  chord  length  required  for  a  given  lift.  The 
center  of  pressure  travel  of  this  type  will  probably  be  large  also  due 
to  the  large  chord  length  and  the  odd  type  of  pressure  distribution  (2) 
with  the  leading  edge  either  inside  or  outside  the  Mach  cone. 

It  is  to  be  noted  that  this  analysis  is  for  the  double  wedge  type 
of  airfoil  cross-sectional  shape  only,  but  the  results  obtained  herein 
will  be  qualitatively  comparaole  for  any  type  of  cross-section. 

Figure  3  shows  the  actual  thickness  ratios  at  which  the  delta  type 
of  airfoil  with  n  *  1.3  and  2.5  will  be  superior  to  the  reverse  arrow 
type  while  Fig.  4  shows  the  same  effect  when  considered  from  the  aero¬ 
dynamic-strength  combination  standpoint.  Limiting  values  of  Kg  at  n  = 

0,  1.3,  1.7,  and  2.5  are  shown  for  purposes  of  comparison.  Fig.  5  shows 
the  maximum  lift-drag  ratios  for  the  range  of  practical  values  of  Eg  for 
values  of  n  *  1.3  and  2.5. 

A  comparison  of  the  maximum  lift-drag  ratios  of  delta  and  reverse 
arrow  wings  is  shown  in  Fig.  7.  The  improved  characteristics  of  the  delta 
type  wing  for  high  thickness  ratios  at  low  Mach  numbers  is  evident  from 
the  curves. 
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In  references  (4)  and  (5),  the  following  expressions  were  derived 
for  lift  and  drag  coefficients: 


Equation  (3)  applies  for  wings  exported  orer  their  entire  hate. 
Kg  will  he  evaluated  helow. 

The  optimum  angle  of  attack  is, 


CCopt  «  28.65  ]/SJc^Z 
and  the  maximum  lift-drag  ratio  becomes 


i  ."W  cl/°l  ~ 

D  max  y  K3  Ai  2  +  * 

Values  of  Cl/Cl^  and  can  he  read  directly  from  Figs.  1  and  2 
respectively.  The  skin  friction  drag  coefficient  is  taken  as  .0053,  as 
in  references  4  and  5.  The  aerodynamic- strength  parameter  £3  is  depend¬ 
ent  on  the  opanwise  pressure  distribution  as  shown  in  the  following  de¬ 
rivation: 


s  »  x  b  x  Z  x  1 

24  5  2  p  o 

where  y/h  is  the  ratio  of  the  di stance  between  the  wing  center  line  and 
the  center  of  pressure,  to  the  wing  semi-span. 

>  By  the  methods  of  reference  (2)*,  it  may  he  shown  that,  for  wings 
whose  leading  edges  are  outside  of  the  Mach  cone, 


♦  This  reference  considers  th®  distribution  over  a  flat  plate  only.  How¬ 
ever  it  i*  assumed  here  that  the  cpanwise  variation  will  he  similar  for 
a  wing  having  a  finite  thickness. 
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where  i  is  found  as  the  limiting  case  of  aquation  ( 28 )  of  that  reference 
to  he 

t  =  I  tan  p  (1  + 

In  the  above  terminology,  is  defined  as  follows: 

Cos/3  =  1 

ir  Tan 

Therefore 

£  =  2  Cos  A  (1  +  2/3  ) 

*  3  7T  Sin  2/3 

The  variation  of  y/b  with/3  as  defined  above  will  be  found  plotted 
in  Fig.  13. 

Now  from  equation  (6) , 

t£  =  y2  =  3  w>^  2  jr 
C2  4  S  b 


Then  X3  =  K4  w  £ 

.  8  S  b 

where  =  G  as  in  references  4  and  5. 

For  the  case  of  the  delta  airfoil  entirely  inside  the  Mach  cone, 
reference  5  indicates  that  the  value  of  y/b  for  large  values  of  n  is 
the  same  as  calculated  above  at  n  =  1.0  (  =0),  namely, 

£  *  JL  =  *4244 
*  sir 

Inasmuch  as  the  slope  of  the  curve  of  Fig.  13  is  0  at  /3  =  0,  it  will 
be  assumed  that  this  value  is  constant  for  any  value  of  n  greater  than 
1.0,  i.e.,  all  cases  of  leading  edge  inside  the  Mach  cone.  Therefore 
the  expression  for  K3  will  be: 


Z3  =  K1  k4  w 
18.85  S 


Since  the  origin.' 1  preparation  of  this  report,  additional  calcul¬ 
ations  for  the  characteristics  of  the  delta  wing  with  the  trailing  edge 
swept  forward  and  backward  h- vc  been  carried  out  by  Dr's.  Puckett  and 
Stow:  rt  (8)  '  n<!  they  have  very  kindly  consented  to  the  inclusion  of 
this  information  in  this  report.  These  data  will  bo  found  in  Figures 
13  through  22.  No  gener;  1  study  of  optimisation  hr.B  been  made  due  to 
the  detailed  investigation  .involved. 
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FORM  DRAG  OF  DELTA  AIRFOILS  WHOSE 
LEADING  EDGES  ARE  INSIDE  THE  MACH  CONE 
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V,  SECTIONAL  CHARACTERISTICS  OF  STANDARD  AIRFOILS 

* 

"by  E.  A.  Borrney 

In  selecting  the  optimum  cross-sectional  shape  of  airfoil  for  use 
in  supersonic  flight,  one  is  tempted  to  choose  the  dou’ole-wedge  type  be¬ 
cause  of  its  low  form  drag  per  unit  thickness.  This  is  a  fallacy  bec'v  -se 
drag  per  unit  stress  is  the  proper  criterion  for  design.  The  increase, 
drag  of  other  types  must  be  balanced  against  their  increased  strength 
(section  modulus)  to  determine  the  optimum  shs  ie. 

With  this  in  mind,  the  following  table  and  curves  have  been  con¬ 
structed  to  show  the  relation  between  drag,  stress,  and  thickness 
ratio  for  various  types  of  airfoils  including  one  designed  by  Dr.  A.  E. 
Puckett  in  which  the  optimum  shape  to  give  minimum  drag  per  unit  stress 
was  the  design  criterion.  It  is  shown  that  the  optimum  shape  is  only 
very  slightly  (0.7  per  cent)  better  than  the  biconvex  (double  circular 
arc)  section.  For  practical  purposes,  therefore,  (inasmuch  as  the 
Puckett  wing  would  be  difficult  to  machine)  the  biconvex  type  of  cross- 
section  is  the  optimum  shape  of  airfoil  for  a  given  allowable  stress  at 
the  root  of  the  wing. 


Shape 

Biconvex 

Optimum 

(Puckptt) 

Double 

VWMtye 
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Double  Wedge 
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Double  Wedge 
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12 

12 

t 

& 

r2 

T2 

r2 

T2 
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24 

12 

I.333T* 
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T2 
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rz 

4 

i 

«!■* 

2a 

wF 
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f 7. 38  9 

24  $ 

18  s 

5  s 

4 

c3 

c3 

c3 

c3 
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&  q  ®  fraction  of  chord  length  having 
wedge  shop*  (each  end) 

O  cn  “  for  infinite  epan  rotio 
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T  =  Vc  (multiply  by  I02  )  T  =  f/c  (multiply  by  ;0"2) 
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